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Significant Findings 

The estimation of microphysical profiles from anvil, convective and updraft regions of an 
oceanic tropical squall were presented. The estimation procedure relied on comparisons 
between brightness temperature calculations and wideband concurrent dual altitude 
brightness temperature observations. Each of the three estimated cloud profiles produced 
calculated brightness temperatures within 10 Kelvin of the observations. It was found that 
a majority of the total iterative adjustments were to the estimated profile's frozen 
hydrometeor characteristics and were required to match the high frequency calculations 
with the observations. The major conclusions of this work include the findings that: (1) 
the microphysical parameters of ice are more detailed and difficult to obtain than are the 
liquid microphysical parameters, (2) cloud resolving models should be validated using 
frequencies >= 150 GHz, and (3) the estimated profiles explain the variability in the 
observations and hence give gross clues as to how the upper cloud is structured. In 
addition, some difficulties matching the 37 GHz observation channels on the DC-8 and 
ER-2 aircraft with the calculations simulated at the two aircraft heights were noted and 
potential causes presented. 



Abstract 


Profiles of the microphysical properties of clouds and raincells are essential in 
many areas of atmospheric research and operational meteorology. In order to enhance 
the understanding of the nonlinear and underconstrained relationships between cloud 
and hydrometeor microphysical profiles and passive microwave brightness temperatures, 
estimations of cloud profiles for an anvil, a convective, and an updraft region of an oceanic 
squall were performed. The estimations relied on comparisons between radiative transfer 
calculations of incrementally estimated microphysical profiles and concurrent dual-altitude 
wideband brightness temperatures from the 22 February 1993 flight during TOGA-COARE. 
The wideband observations HO — 270 CW*) nm , , 
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reaching up to 20 km. The low frequencies enhance the rain and cloud water profiles, 
while the high frequencies are required to detail the higher altitude ice microphysics. A 
microphysical profile was estimated for each of the three regions of the storm. Each of 
hree estimated profiles produced calculated brightness temperatures within -~10 K 
of the observations. A majority of the total iterative adjustments were to the estimated 
profile's frozen hydrometeor characteristics and were necessary to match the high frequency 
calculations with the observations. This indicates a need to validate cloud resolving models 
using high frequencies. Some difficulties matching the 37 GHz observation channels on 
the DC-8 and ER-2 aircrafts with the calculations simulated at the two aircraft heights 
(~ 11 km and 20 km, respectively) were noted and potential causes presented. 
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Popular Summary 


Profiles of the microphysical properties of clouds and raincells are essential in many 
areas of atmospheric research and operational meteorology. In order to enhance the 
understanding of the nonlinear and underconstrained relationships between cloud and 
hydrometeor microphysical profiles and passive microwave brightness temperatures, 
estimations of cloud profiles for an anvil, a convective, and an updraft region of an 
oceanic squall were performed. The estimations relied on comparisons between radiative 
transfer calculations of incrementally estimated microphysical profiles and concurrent 
dual-altitude wideband (10—220 GHz) brightness temperatures observed during TOGA- 
COARE The low frequencies enhance the rain and cloud water profiles, while the high 
frequencies are required to detail the higher altitude ice microphysics. Each of the three 
estimated cloud profiles produced calculated brightness temperatures within 10 Kelvin of 
the observations. 
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1. Introduction 

An understanding of the relationships between passive microwave brightness 
temperatures (T B ) and the microphysical characteristics of hydrometeors has long been a 
goal of atmospheric scientists (Wilheit et al. 1982; Smith et al. 1992; Wang et al. 1998). 
Knowledge of these relationships is essential for many reasons including: developing cloud 
resolving models (e.g., Tao and Simpson 1993), improving forward radiative transfer 
calculations (e.g., Smith et al. 1992), computing vertical diabatic heating profiles (e.g., Tao 
et al. 1993), and in precipitation parameter retrieval algorithms (e.g., Marzano et al. 1999; 
Smith et al. 1992). In order to gain information about these relationships many theoretical 
(e.g., Adler et al. 1990: Smith et al. 1992) and experimental fe.ff.. Hevmsfield et al 199fv 
Simpson et al. 1996; Liu and Curry 1998) studies have been performed. In general, these 
studies compare T B observations to T B calculations using estimated microphysical profiles 
or use empirical knowledge gained from experimental observations. The results from these 
studies have clarified the relationships between the size distributions of liquid (rain and 
cloud water) hydrometeors and the upwelling brightness temperatures over oceans. While 
the liquid hydrometeor relationships are fairly well understood for low frequencies (Wilheit 
1986), with only minor differences for varying cloud types, the relationships between the 

higher altitude frozen hydrometeors in convective clouds and the brightness imagery is less 
clear. 


The relationships between frozen hydrometeors and the T B are more difficult to 
grasp for several reasons, including inherent nonlinearities, the enormous variety of frozen 
particle shapes, sizes, and air-ice-water ratios (LaChapelle 1969), the sparsity of measured 
and/or "truC’ ice profiles along with coincident, colocated T B observations and the low 
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sensitivity of the commonly-used lower frequencies to the frozen hydrometeors. Many 
theoretical studies have begun to address the nonlinear relationships between the T B and 
the shapes and sizes of frozen hydrometeors (e.g., Evans and Stephens 1995). In addition, 
several experimental field campaigns have attempted to coincidently measure both T B s 
and cloud microphysical characteristics (e.g., Webster and Lucas 1992; Simpson et al. 1996; 
McFarquhar and Heymfield 1996) with limited success. Many of these experimental studies 
use low to mid frequencies (< 90 GHz) to investigate and estimate microphysical profiles. 
However as shown in this paper, many estimated profiles with different ice microphysics 
can match the < 90 GHz frequency observations. 

By using both low and high frequency observations, a more accurate relationship 
between the T B and the liquid and frozen hydrometeors can be found. The low frequencies 
can be used to estimate the low altitude liquid hydrometeors, while the high frequencies are 
necessary in enhancing the higher altitude ice microphysical details (Evans and Stephens 
1995). This paper will describe how coincident wideband (10—220 GHz) passive microwave 
T b observations from radiometers on both the DC-S aircraft (~ 11 km altitude) and 
the ER-2 aircraft (~ 20 km altitude) obtained during TOGA-COARE (Tropical Ocean 
Global Atmosphere/Coupled Ocean Atmosphere Response Experiment) were utilized to 
investigate the underlying microphysical cloud profiles. A nadir-viewed, dual-altitude 
comparison between TOGA-COARE brightness temperature observations and computed 
brightness temperatures was the basis for the microphvsical profile estimation procedure. 
Microphysical profiles used in the calculations were initialized with profiles obtained from 
the Goddard Cumulus Ensemble (Tao and Simpson 1993) cloud parameterization model. 

Ice densities and drop size distributions were augmented with densities obtained during 
experimental observations. By systematically varying the microphysical cloud parameters 
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(e.g., particle size distributions and the hydrometeor density) and requiring that the 
calculations remain consistent with the wideband TOGA-COARE observations, the cloud 
conditions were estimated for three regions of an oceanic convective storm. 

Microphysical profiles for an anvil, convective, and convective updraft region 
of a storm observed on 22 February 1993 during TOGA-COARE are estimated. While 
knowledge of the estimated profiles for the three case specific regions reported herein 
is interesting, this work identifies the advantages of wideband and in particular, high 
frequency (> 90 GHz) observations. High frequency observations can clarify high altitude 
(> 8 km) microphysics. Improved knowledge of ice microphysics can then be used to 
validate and enhance cloud resolving models. Cloud resolving models must arronnt for Wh 
low and high altitude microphysics in order to accurately represent cloud systems. This 
work compliments current cloud model verification efforts which focus on low frequency 
observations and hence low altitude microphysics. Another advantage of a detailed 
wideband comparison is to validate instrument calibration. 

This paper begins by describing the observational data from TOGA-COARE. In 
particular, the instruments and the three regions of interest for a flight on 22 February 1993 
are detailed. The T B calculations and baseline vertical microphysical profiles are described. 
Then the procedure used to compare the observations with the calculations in estimating 

the microphysical profiles for the three regions is outlined. The last section summarizes the 
conclusions. 
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2. Observations 

The brightness temperature observations of interest were obtained early in 1993 
during TOGA-COARE. Among the sensors that measured the response to atmospheric 
hydrometeors in the microwave region of the electromagnetic spectrum during TOGA- 
COARE were the AMPR (Advanced Microwave Precipitation Radiometer) (Spencer et al. 
1994), and MIR (Millimeter-wave Imaging Radiometer) (Racette et al. 1996) on board the 
ER “ aircra ft> and the ARMAR (Airborne Rain Mapping Radar) (Durden et al. 1994), and 
AMMR (Airborne Multi-channel Microwave Radiometer) (Wang et al. 1994), and AMMS 
(Airborne Microwave Moisture Sounder) (Wang and Chang 1990) on board the DC-8. The 
MIR observation frequencies were 89, 150, 183.3T1. 183.3T3. 183.3T7. and 220 GHz while 
the AMPR observed at 10.7. 19.35. 37, and 85 GHz. On the DC-S aircraft the AMMR 
operated at 18 . 1 , 21, 37, and 92 GHz, and the frequencies for the AMMS were 92, 183. 3±2, 
183.3±o, and 183. 3±9. The footprint sizes ranged from 870 m to 2.8 km for the AMMS 
183.3 GHz channels and the AMPR 10.7 GHz channel, respectively. The ARMAR was a 
single-channel 13.8 GHz radar. 

There were a few coordinated flights where the DC-S under-flew the ER-2 and 
provided a rich data set of observed brightness temperatures at 2 altitudes (~ 11 km DC-8 
and - 20 km ER-2). One such coordinated data flight occurred on 22 February 1993. 

Here the DC-8 flew through an oceanic convective storm and the ER-2 flew over the same 
latitude and longitude (fairly straight flight track from 9.503°S, 15S.577°E to 9.503*S, 
lo9.0S2 c 'E) just a little over 2 minutes later. This well studied event (e.g., Wang et al. 1998; 
Jorgensen et al. 1995; Trier et al. 1996; Marzano et al. 1999) was near Guadalcanal Island 
in the Western Pacific with a north-south oriented squall-line. Figure 1 shows the radar 
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reflectivities from the ARMAR for this flight track. The convective system was slanted 
from right to left with increasing altitude on the left hand side. The freezing level was near 
4 km altitude, however high reflectivity values exist above the freezing level. In the storm 
core, some radar reflectivities exceed 45 dB or are severely attenuated (green-colored area). 

Figure 2 shows the low and mid-frequency (< 90 GHz) nadir-viewed T B values 
for each aircraft, while Fig. 3 shows the high frequency T B values. These figures reinforce 
the qualitative information obtained from the radar reflectivity plot in Fig. 1. Figure 2 
shows evidence of a scattering signature in the 37 GHz channels from about 21:32:40 (UTC 
DC-8). There is also evidence of scattering at 19 and 21 GHz near 21:34:40 (UTC DC-8) 
where there is a slight depression in these T B values. For Fie. 3. the left hand tilt and 
anvil region (arrow A) of the storm becomes noticeable in the differences between the DC-8 
and ER-2 T B values where there is a 100 K difference due to ice scattering at the upper 
altitudes. An area of uniform convection is identified with the C arrow. Note that the 
10./ GHz channel only reaches a maximum near the convective updraft region (arrow U) 
where there is significant liquid water as indicated by the attenuation from rain in the radar 
reflectivities (green-colored). 

The first region of interest is the anvil region (21:33:14 — 21:33:55 ER-2 UTC 
21:31:08 — 21:31:49 DC-S UTC, arrow A). This region is characterized by a theorized 
high ice cloud above the ARMAR field of view. Comparing the DC-8 and ER-2 183 GHz 
channels in Fig. 3 (arrow A) shows that the ER-2 183 GHz channels are about 100 Iv lower 
than the DC-S 183 GHz channels. These significantly lower T B values along with 90 GHz 
channels with similar responses (Fig. 3) indicate that frozen hydrometeors at altitudes 
> 10 km (an anvil cloud) are effective scatters at frequencies greater than 90 GHz. 



7 


The second region of interest is a convective region (21:35:37 — 21:35:48 ER-2 
UTC. 21:33:07 — 21:33:18 DC-8 UTC). This region (arrow C in Figs. 1, 2, 3) is characterized 
by 45+ dB radar reflectivity values from about 2 to 4 km and fairly significant radar 
reflectivity returns above the freezing level. The convection is fairly uniform in this region 
and the T B values show that the 183.3 GHz channels are about 60 I\ cooler for the MIR 
than for the AMMS on the DC-8 indicating the existence of ice scatters above the DC-8 
aircraft. The warm low frequency T B values (Fig. 2) are due to the presence of liquid (rain 
and cloud water) absorbers below the freezing level. 

The third region (21:37:00 — 21:37:06 ER-2 UTC, 21:34:24 — 21:34:32 DC-8 
UTC) is characterized by strong updrafts and high radar reflectivities. This reprion (arrow 
U in Figs. 1, 2, 3) is not uniform and shows significant vertical and horizontal variability. 
The radar reflectivities (fig. 1) show significant attenuation at low altitudes caused by 
heavy rainfall. In Fig. 2 the low frequencies are warm indicating the presence of rain, while 
the high frequency T B values (Fig. 3) are very low suggesting large amounts of ice. 

Other experimental data used in this research occurred during the Central 
Equatorial Pacific Experiment (CEPEX) in March and April of 1993 (McFarquhar and 
Heymfield 1996). CEPEX sampled frozen hydrometeor number densities and drop size 
distributions in the same Australian warm pool region as during the TOGA-COARE 
experiment observed. Liu and Curry analyzed this data and presented profiles of cloud ice 
average size, average density, and number density 1 as a function of temperature (Liu and 
Curry 1998). As will be shown in Section 4, the CEPEX data is used to enhance the ice 

'Note that the units of A' 0 in Fig. 2 of (Liu and Curry 1998) are incorrectly specified. 
The correct units are m -3 ^m _1 (Liu 1997). 



8 


microphysics in the estimated profiles. 

3. Brightness Temperature Calculations 

A radiative transfer model is required to transform the (iteratively) estimated 
microphysical atmospheric profile information into upwelling passive microwave brightness 
temperatures. The planar-stratified, scattering-based radiative transfer (RT) model used 
herein was originally developed by Gasiewski and Staelin (1990) and later modified to 
allow for five (or more) hydrometeor types (e.g., suspended cloud water, rain, suspended 
cloud ice, snow, graupel) (Skofronick-Jackson and Gasiewski 1995). The radiative transfer 
model requires as input vertical profiles of temperature, height, relative humidity, and five 
hydrometeor parameters (per hydrometeor type): type, density (gm' 3 ). particle average 
radius ia 0 , mm), particle intercept ( A 0 , mm^m -3 ), and a flag identifying if the particle 
has a monodisperse or polydisperse size distribution. Flexibility exists in that the user can 
select observation height, viewing angle, frequency and polarization. In this work a nadir 
viewing angle is assumed and the observation heights and frequencies match those of the 
ER-2 and DC-8 aircraft instrumentation. Since this work uses nadir Jgs observed with a 
narrow beam width, the plane parallel radiative transfer model is acceptable. In addition, 
a calm ocean surface with a 27°C temperature and with wind speeds < 6 ms -1 is assumed. 
The ocean is obscured for all but the lowest frequencies for all three storm scenarios, thus a 
calm ocean model is acceptable. 

In computing the absorption and scattering coefficients for the individual 
hydrometeors, decaying exponential polydispersive size distributions are assumed (unless 
the drop size distribution dictates a monodispersive size distribution). The RT code uses 
either Rayleigh or Mie theory based on the electrical size of the hydrometeor. When the 
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electrical size is small (|i>y| < 0.5, where n is the complex refractive index and x = 2e 
(U laby et al. 1981)) the easier-to-implement Rayleigh theory is used. Typically dielectric 
mixing theory (Borhen and Battan 1980) is used to determine the effective dielectric 
constant (permittivity, n 2 j of the heterogeneous hydrometeors (e.g., snow and graupel). 
The fractions of entrained ice, air, and water for these particles must be specified along 
with their drop size distributions. The validity of this (and most) dielectric mixing theories 
is void when the size of the inclusion is greater than or equal to the wavelength (Borhen 
and Battan 1980; Sihvola 1989). For a typical mid-latitude oceanic convective storm, 
validity is ensured to at least 100 GHz (Skofronick-Jackson and Gasiewski 1995). Since this 
research includes frequencies up to 220 GHz, all frozen hydrometeors have been restricted 
to homogeneous solid ice so that dielectric mixing theory is not needed. In addition, a 
spherical shape is assumed for all hydrometeors. 

The estimated microphysical profiles are initially obtained from the Goddard 
Cumulus Ensemble (GCE) three-dimensional simulated cloud model (Tao and Simpson 
1993). The GCE simulation used in this research was initialized using nearby composite 
aircraft and radiosonde sounding obtained a few hours earlier than the DC-S/ER-2 
observations on 22 February 1993. The GCE model produces 4-dimensional distributions 
of hydrometeors that have been used as input to passive microwave radiative transfer and 
space-borne precipitation radar models by researchers for improving their precipitation 
retrievals (e.g., Prasad et al. 1995: Yeh et al. 1995; Simpson et al. 1996). For the 22 
February 1993 GCE cloud model, the assumed size distributions, or number of particles 
with diameter D, follow a decaying inverse exponential drop size distribution: 

h h{D) = Nh 0 e 2oft o m -3 mm -1 


( 1 ) 
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where (Rutledge and Hobbs 1984) 


ah o = 2.0 


M h 

JPhNh, 


1 0.25 


= (r) mm. 


■0 J 


( 2 ) 


In the above. M h is the partial density (gm“ 3 ) of hydrometeor h in a cubic meter of air, p h 
is the average intrinsic density (gm -3 ) of the hydrometeor, N ho (m -3 mm -1 ) is the intercept 


multiplier and (r) is the average hydrometeor radius. The partial density M h is provided 
by the GCE or CEPEX microphysical cloud data. Here h is used to distinguish among 
the various classes of hydrometeors (e.g., p w , p r , p { , and p g denote the intrinsic density 
of cloud water, rain, ice, snow, and graupel, respectively). The rain, snow, and graupel 
hydrometeors have the exponential size distributions of Rutledge and Hobbs (1984) given 


m Eqno. 1 mid 2 with parameters N 0h = 8000.0, 4000. u. 4uuu.u m 3 mm ! respectively and 
the specific density ph = 1.0 x 10 9 , 1.0x10 s , 4.0x10 s gm -3 , respectively. For the referenced 
GCId case, the ice-air-water ratios for cloud water, rain, cloud ice, snow, and graupel are 
0-0-10026, 0-0-100%, 100-0-0%, 10-90-0%, and 40-60-0% respectively (Tao and Simpson 
1993). Note that there is no specified melting layer in the GCE model. The suspended 
cloud water and ice particles are small enough to advect with the airflow and have a fixed 
mean diameter of (r w ) = 0.003 mm and (r.) = 0.02 mm. Here, the parameters N w0 and N l0 
(and thus the number of particles) vary to account for the differing densities M w and M t . 


The GCE cloud model is typically used to compare with radar and brightness 
temperature observations (e.g., Adler et al. 1991; Marzano et al. 1999). The emphasis 
in these comparisons, however, is based on a general gross scale and on the inherently 
low altitude rain structures. There were more degrees of freedom when comparisons were 
made in the regions of mixed-phase and ice clouds at higher altitudes. Since the scattering 
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coefficient in the Rayleigh region goes as (Borhen and Huffman 1983) 

~ JVo/ 4 r 6 (3 

(where / is the frequency), the size (r) of the particle dominates the scattering term. 
Thus small changes in r can lead to significant changes in the brightness temperature. In 
the Mie region, the relationships are less clear, however, particle size is still extremely 
important. In general, the GCE modeled particles do not reproduce the high frequency 
brightness temperatures observed by the airborne radiometers for the three regions of 
interest from the 22 February 1993 TOGA-COARE flights. This suggests the need of 
adjusting the model-produced cloud microphysics and possibly modifying the parameters 
in the initialization of the OCR modH. 

In order to adjust the GCE produced microphysical profiles and to provide 
additional hydrometeor types, two other drop size distributions were included in the 
radiative transfer code. The CEPEX data, as analyzed by Liu and Curry (1998), 
was incorporated into the code. However because the CEPEX data is derived from 
measurements (McFarquhar and Heymfield 1996) there is not a realistic method for varying 
the density, or drop size in order to vary the microphysical profile and thus the calculated 
T b . In order to provide more flexibility in the frozen hydrometeor drop size distributions, 
the Sekhon and Srivastava (1970) drop size distribution was also included. Generally the 

CEPEX density was used in the Sekhon-Srivastava (SS) drop size distribution equations 
which follow Eqs. 1 and 2 with 

a 0 = 0.422 M 0 - 52 mm, ( 4 ) 

and 


Ao — 636 .' AM 109 mm l m 3 . 


( 5 ) 
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Because the SS drop size distribution is of a functional form, changes in M, a 0 ,or N 0 can be 
applied and propagated through to the absorption and scattering coefficient calculations. 
Changes in M or ao were often necessary because the CEPEX data was an average of 
the measured microphysics and the cases studied in this work were at the extreme ends 
of convective and anvil cloud structures. A comparison between the CEPEX drop size 
distribution and the SS drop size distribution shows that for a fixed density M, the SS 
average radius is smaller and N 0 is larger than for the CEPEX case. 

4. Estimations 

A wideband, nadir-viewed, dual-altitude comparison between TOGA-COARE 
brightness temperature [Tb) observations and computed brightness temperatures was the 
basis for the estimation procedure. A feedback procedure in the comparison process is 
used to iteratively improve the estimated microphysical profiles. The feedback procedure, 
driven by an understanding of hydrometeor — Tb relationships, adjusts the microphysical 
parameterizations of hydrometeor density, drop size distribution, and the ice- air- water 
fractions of the frozen hydrometeors. In order to obtain the best retrieved profile, changes 
in the hydrometeor parameterizations were not constrained to conserve the physical 
parameters (e.g., mass) of the initializing profile. By varying the microphysical cloud 
parameters (e.g., particle size distributions and the hydrometeor density) and requiring 
that the comparisons remain consistent with the wideband TOGA-COARE brightness 
temperatuie observations, the cloud conditions between the two aircraft were estimated. 
Concurrent radar ieflecti\ ities from the ARMAR on the DC-S are also used in the analysis 
to provide a crude profile of rain densities when the GCE data was inadequate. 


The feedback/iterative process is stopped when at least fifteen out of sixteen 
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observed microwave channels of the ER-2 and DC-8 are within ±10 K of the calculated 
Tb values. While not adopted in this work, the ±10 K requirement could be lessened 
for channels with larger footprint sizes due to the microphysical information lost during 
the observation process. Furthermore, the RMS error of the calculated minus observed 
differences was minimized. While the estimated profiles do match the T B observations fairly 
well, there is no guarantee that these profiles are the global minimum in the solution space. 
While the retrieved profiles may not be unique, variations of the profiles would only affect 
the RMS error. The conclusions are not expected to change. Furthermore the estimated 
profiles do explain the variability in the observations and hence give gross clues as to how 
the upper cloud is structured. 

(a) Anvil Region 

For the anvil case, the 12Sxl2S profiles of the most mature time frame (360 
minutes) of the three-dimensional GCE microphysical cloud were sorted to determine which 
profiles contained large amounts of frozen hydrometeors at higher altitudes (> 11 km). 
Eight favorable vertical profiles were selected and extracted for use as baseline microphysical 
profiles. Each of the eight profiles were non-precipitating with frozen hydrometeor densities 
exceeding 2 gm 3 above 11 km. (These were the maximum upper altitude ice densities in 
the available GCE profiles.) Using the drop size distributions suggested in Tao and Simpson 
(1993) and described in Section 3 the brightness temperatures were computed for the eight 
piofiles. The brightness temperatures were averaged together to compare with averaged 
brightness temperatures from the specified anvil region UTC times of the ER-2 and DC-8 
(arrow A in Figs. 1, 2, 3). 


The resultant averaged brightness temperatures are provided in Table 1. In 
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comparing these values with the observations, all calculated values for frequencies 
> 150 GHz are too warm by up to SO K with respect to the ER-2 observations at those 
frequencies. On the other hand, the high frequency DC-8 altitude calculations are cooler 
than the DC-8 observations. In fact, the calculated T B values are too cool by 10-50 I< for 
all the DC-S observed frequencies. These comparisons indicate several items about the 
calculations: (1) there is not enough small-sized ice at the higher (> 11 km) altitudes, 

(2) that the ice-air mixtures in the snow and graupel do not produce enough scattering, 

(3) that dielectric mixing theories used for the GCE snow and graupel may be breaking 

down at the highest frequencies, or (4) there may be too much ice, or incorrectly-sized ice 
at the lower altitudes. 

In attempting to correct for the inconsistencies described above, the frozen 
hvdrometeor components of the GCE profile were replaced by averaged ice densities and 
drop size distributions as derived by Liu and Curry (1998) using the CEPEX data. The 
selected GCE profiles were modified such that all the cloud ice, snow, and graupel were 
replaced with the CEPEX densities and drop size distributions whenever T < 273.15 K, 
(i.e., at heights >~ 5 km). The resultant T B values are provided in Table 1. This 
table shows that while the ER-2 simulated calculations are closer to the observations, 
the DC-8 calculations are in worse agreement with the observations. These unfavorable 
comparisons indicate that the CEPEX densities and drop size distributions will require 
further adjustment prior to obtaining a favorable comparison. 

After evaluating many ice density and drop size distributions through a manual 
feedback analysis, the best microphysical parameterization used a hybrid profile (Fig. 4) 
combining the GCE profile, CEPEX ice. and variations of the Sekhon and Srivastava (1970) 
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drop size distribution. The calm ocean surface temperature was assumed to be 301 K and 
the relative humidity varied from 73% to 100%. In the hybrid profile, the GCE cloud water 
densities were replaced with the GCE rain water densities (which were small to start with) 
and the rain densities zeroed. Because the cloud water is monodisperse with a fixed drop 
diameter of 0.015 mm, less absorptive warming occurs and the 18 - 89 GHz calculations 
match better with the observations. With the zeroed rain densities, the calculated 10.7 GHz 
channel was still too cool, thus the cloud water densities below 250 m were multiplied by 
4.0 to add absorptive warming to the 10.7 GHz channel. The frozen hydrometeors started 
at an altitude of 8 km (-15°C) where the CEPEX ice densities were used with the SS drop 
size distribution up to an altitude of 16 km. Starting the ice at a higher altitude than 
' llC flfcciill o lcVeI (~ S kni) is more appropriate given the low return in the low altitude 
radar reflectivities of Fig. 1 and the too cool DC-S T B calculations of Table 1. From 12 
to 14 km, the CEPEX ice densities using the CEPEX drop size distribution is added to 

increase the scattering component thus reducing the high-altitude high-frequency brightness 
temperatures. 

The resultant brightness temperatures for the modified drop size distributions are 
provided in Table 1. In comparing these with the observations on 22 February 1993 from 
21:33:14 — 21:33:49 UTC (ER-2), the results are favorable. All T B values are within ±10 I< 
of the observations. The average error (absolute value of the differences) of the hybrid 
profile is 4.5 K, while the RMS error of the differences is 1.3 K. 

(b) Convective Region 

For the convective case, the GCE data was used to provide a baseline profile of 
height, temperature, relative humidity, and cloud water density. The GCE data was not 
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used to provide profiles of frozen hydrometeor densities due to dielectric mixing issues 
associated with the GCE inhomogeneous snow and graupel hydrometeors. In addition, 
instead of sorting through the 128x128 GCE pixels for a suitable rain profile, a crude 
baseline rain water density was determined using the reflectivity profiles from the ARMAR 
(arrow C in Fig. 1). Initially, no frozen hydrometeors were present — they were added 
during the feedback procedure. The relative humidity of the GCE profile was reset to 100% 
whenever hydrometeors were present in the profile. 

In order to obtain matching computed T B values, several modifications to the 
baseline profile were necessary. First the rain water required adjustment so that the low 
frequency calculations better matched the observations. The adjustment increased the rain 
density by 10% below the freezing level 5 km) and decreased the rain density by 20% 
abo\e the freezing level (Fig. 5). This adjustment brought the calculations for < 19.35 GHz 
to within 8 k of the T B observations. Since no frozen hydrometeors were initially present, 
the next step was to add frozen hvdrometeors to the lower levels (0°C to ~ 11 km) in 
order to get the high-frequency low-altitude calculations to match the DC-S high frequency- 
observations. To achieve a better representation, the CEPEX ice densities with a modified 
SS drop size distribution were added above 6 km. The SS drop size distribution was 
modified by multiplying the average radius by 1.475 in order to increase the scattering 
component. In performing the drop size distribution modification, the number density was 

kept at the original SS value given in Eq. 5. Thus the actual density was modified from the 
CEPEX density by: 

Mnew = McEPExC 4 (6) 

where C is the constant multiplied by « 0 (e.g., 1.475). Above the height of the DC-8 frozen 
hydrometeors with a smaller average radius were required to obtain a good match with 
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the ER-2 high frequency observations. Thus above 11.3 km the average SS drop radius as 
generated from the CEPEX density was divided by 1.5. These small-sized hydrometeors 
were added to a height of 16 km. In addition, to increase the scattering component, CEPEX 
densities using the CEPEX drop size distribution was included from 11.3 to 15 km (Fig. 5). 
A final adjustment to change the relative humidity to 100% whenever hydrometeors were 
present further improved the water vapor channels. 

These modifications to the original profile brought 15 of the 16 T B calculations to 
within 10 K (Table 2) for all observations (arrow C in Figs. 1, 2, 3). Potential causes for 
the mismatched ER-2 37 GHz channel are provided in Appendix A. The average absolute 

value difference between the observations and calculations is 4.273 K while the RMS error 
is 5.5S7 K. 

(c) Convective Updraft Region 

1 he radar reflectivity profiles associated with the convective updraft region (Fig. 1 
anow l) indicate strong updrafts lifting rain and ice to altitudes greater than 9 km. There 
is no clearly defined melting layer, however low altitude attenuation implies heavy rainfall. 
Thus in defining the baseline profile, the radar reflectivity profile was used to obtain a 
profile of rain densities. The temperature, height, and cloud water profiles were extracted 
from a single GCE pixel with convective characteristics. The relative humidity was set 
equal to 100% whenever hydrometeors were present. 

Since our simple radar reflectivity to rain density conversion did not account 
for attenuation, additional ram water was required to warm the low frequency channels. 
Because of the updraft nature of this region the rain water extended above the freezing 
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level to approximately 9 km (Fig. 6). This brought the low frequency channels to within 
7 K of the radiometer observations on both the ER-2 and DC-S aircraft (Table 3). 

The next goal was to obtain a match to the DC-8 observations sensitive to the 
frozen hydrometeors. Initially the CEPEX ice densities and drop size distributions were 
used for the frozen hydrometeors. The brightness temperatures resulting from the addition 
of the CEPEX ice were too warm by an average of 19 K at 37 and 92 GHz and too cool by 
an average of 15 K for the 183.3 GHz channels. These unmatched brightness temperatures 
imply that there is not enough scattering at the altitudes where the 37 and 92 GHz channels 
are sensitive, while scattering is too intense for the 183.3 GHz channels. Thus several 
modifications were made to the microphysical profile. 

First to correct for the warm 37 and 92 GHz calculations, we note that Gasiewski 
(1993) shows that as hydrometeor density increases (over a uniform rain cell) the 
temperature weighting functions peak at increasing altitudes. Thus for this hydrometeor 
rich convective updraft region we increased the scattering component at altitudes high for 
the DC-8, e.g., from 8.25 11 km. In this range, the CEPEX density was multiplied by 

11.75 and the SS drop size distribution was used (See Fig. 6). The ice below 8.25 km was 
left with the CEPEX drop size distribution. This change brought all the calculated DC-8 
altitude brightness temperatures to within 9 K of their corresponding observations. The 
low altitude 183.3 GHz Tb calculations are warmed because the heavy ice layer between 
8.25 and 11 km provides opacity such that these channels do not probe as deep into the 
frozen hydrometeor layers. Furthermore, the relative humidity of the profile was increased 
to 100% whenever hydrometeors were present. 


The final step in the estimation procedure was to cool the high altitude (ER-2) 
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high frequency Tg calculations by adding ice to the upper altitudes. The microphysical 
profile below 11 km was retained in order to preserve the DC-8 Tg calculations. Above 
11 km the CEPEX density using the CEPEX drop size distribution was used to the top 
of the CEPEX cloud top (at 15 km). CEPEX ice alone did not provide enough cooling. 
Therefore a second ice layer was added where the CEPEX density profile was multiplied 
by 0.2o and used in the SS drop size distribution from 11-16 km (See Fig. 6). These 
modifications to the upper altitude ice brought all but the ER-2 37 GHz channel to within 

9 k of the observations (See Table 3). Several causes for difficulties at 37 GHz are noted in 
Appendix A. 

5. Conclusions 

The estimation of microphysical profiles from anvil, convective and updraft 
regions of an oceanic tropical squall have been presented. The estimation procedure 
relied on comparisons between Tg calculations and wideband concurrent dual altitude 
brightness temperature observations. There are several major conclusions of this work: 

(1) the microphysical parameters of ice are more detailed and difficult to obtain than 
are the liquid microphysical parameters, (2) cloud resolving models should be validated 
using frequencies > 150 GHz, and (3) the estimated profiles explain the variability in the 
observations and hence give gross clues as to how the upper cloud is structured. To support 
these conclusions, the calculations initially assumed hydrometeor microphysics as given 
in a widely-accepted cloud resolving model. In all three regions of the storm, the cloud 
model microphysics required substantial modification in order to obtain an acceptable 
match with the observations. It was found that low frequency (< 37 GHz) matching 
was easily accomplished with minor variations in the liquid hydrometeor densities and 
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drop size distributions. However, for frequencies > 150 GHz the brightness temperature 
observations required detailed frozen hydrometeor microphysics and approximately 80% 
of the feedback iterations were devoted to achieving T B calculations that matched the 
high frequency observations. We found that the more complex the vertical profile, the 
more difficult it was to minimize the RMS error and to obtain less than 10 K calculation 
minus observation differences. Due to the significant effort needed to obtain a match to 
the higher frequencies, and the fact that multiple ice microphysical profiles generated T B 's 
that matched observations < 90 GHz, cloud resolving models should be validated with 
observations at frequencies > 150 GHz. 

This work also showed that the three storm scenarios required very different 
ice microphysics in order to produce T B calculations matching the TOGA-COARE 
observations. By comparing Figs. 4-6, it can be seen that the maximum ice density 
(0.07 gm- 3 ) and average drop size (0.11 mm) in the anvil profile are much smaller than for 
the convective and updraft cases. For the convective and updraft cases, the maximum ice 
density in the profile is nearly equivalent at 0.51 and 0.48 gm' 3 , respectively. However, 
the average ice particle sizes for the updraft profile are double those of the convective case. 
Since the radiative transfer model extends the polydispersive exponential drop sizes to 1 1 
times the average size (Gasiewski 1993), the maximum drop radii are 1.2, 2.1, and 4.0 mm 
for the anvil, convective, and updraft regions, respectively. 


While this paper presents an improvement in estimating high altitude ice 
microphysics by using high frequency information, there are at least two approaches that 
can be used in the next generation of this estimation. First, quantitatively incorporating 
radar reflectivity observations from the ARMAR and or EDOP would further improve the 
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microphysical profile retrievals. In fact, Meneghini et al. (1997) and Marzano et al. (1999) 

have already started this type of comparison. In addition to make the microphysical profiles 

more realistic, a layer of inhomogeneous ice particles (a melting layer) should be added. 

However such a melting layer must be mindful of the frequency— inclusion size limitations 

associated with the dielectric mixing theories. Thus the size of the inclusions must be very 

small and/or the melting layer placed at an altitude where the high frequencies do not 
probe. 


Appendix A: Difficulties at 37 GHz 


With respect to the difficulties in matching the 37 GHz channels, we note that 
tor both convective and convective updraft regions the observed T B for the ER-2 is warmer 
than the DC-3 T B values at 3/ GHz. Warmer ER-2 T B values are not seen in the anvil 


region. There are a few factors that may contribute to these difficulties. These factors 
are of three types: (1) temporal flight displacement issues, (2) modeling problems, and 
(3) estimation algorithm effects. These factors are discussed in this appendix, however, 

all of these causes are difficult to verify and suggest the complexity of studying the storm 
structure in detail. 


same 


First the 2-3 minute temporal displacement in the observations over the 
locations between the ER-2 and DC-S aircraft could present problems in these comparisons 
because of the fast moving and changing convective core structure. A close examination of 
the AMPR image in the neighborhood of the convective updraft region shows that a cold 
cell with 3t GHz T b values of ~ 22-5 K appears about 2 km upstream and a couple of pixels 
from nadir. This could be the same cell observed by AMMR on board the DC-S aircraft; a 
positive identification is not possible because AMMR is not an imager. On the other hand, 
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the anvil region is expansive and characterized by a gradual and slow change in T B , and is 
therefore unlikely to be subject to the same difficulty. 

Next, reproducing the effects of the microphysical parameterization model in the 
radiative transfer model can be difficult at 37 GHz. The 37 GHz channel is much more 
sensitive to the ratio of frozen to liquid hydrometeors because it responds to both the liquid 
and frozen hydrometeors fairly equally (whereas, 19 GHz responds mostly to liquid and 
only slightly to frozen hydrometeors and 89 GHz responds more to ice than liquid). This 
may mean that our retrieved profile still does not have the proper ratio of frozen to liquid 
hydrometeors (above the DC-8 altitude as this is where the 37 GHz calculations are too 
cool) or that the ratio of frozen versus liquid hydrometeors has changed in the 2-3 minute 
delay between the DC-S and ER-2 flights. In addition, all frozen hydrometeors have been 
parameterized as solid ice spheres. In reality, it is likely that in the convective and updraft 
regions the ice will be an inhomogeneous mixture of ice, air, and water. Adding a melting 
layer with inhomogeneous ice would warm the low and mid-frequency T B calculations due 
to the increased absorption component. A concern with adding a melting layer is that 
the dielectric mixing theories needed to compute the absorption and scattering coefficient 
require the size of inclusion to be much smaller than the wavelength. This will not be 
the case for the higher frequencies. In fact, initial tests show almost a 1 K change in T B 

for > 150 GHz and more change for lower frequencies when the ice between 0 and -7° C 
contains 10% liquid water. 

A third factor is that the estimated microphysical profile may only be a local 
minimum in the solution space and other microphysical profiles might be more appropriate. 
Finally, the AM MR measurement accuracy is about ±4 Iv (Wang et al. 1994). It is quite 
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possible that the AMPR calibration accuracy is of the same order of magnitude and thus a 
measurement discrepancy of ~ 5 — 6 K between these instruments is not totally unexpected. 
If there are measurement discrepancies for the 37 GHz channels, the estimation algorithm 
should reduce the ±10 K requirement for 37 GHz. 
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Figure 1: The ARMAR (13.8 GHz) radar reflectivity profiles for 22 February 1993. The 

ARMAR was on the DC-8 and measured reflectivities to 9 km. The UTC times are of the 
DC-8 flight. 


Figure 2: The low frequency brightness temperature values for 22 February 1993 as observed 
from the (a) ER-2 and (b) DC-8. The UTC times are of the ER-2 (a) and DC-8 (b) flights. 

Figure 3: The high frequency brightness temperature values for 22 February 1993 as 

observed from the (a) ER-2 and (b) DC-8. The UTC times are of the ER-2 (a) and DC-8 
(b) flights. 
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Figure o: The estimated profile for the convective region: (a) densities, (b) average radii. 
Figure 6: The estimated profile for the convective updraft region: (a) densities, (b) average 
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Table 1: Anvil Tg Comparisons 
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Table 2: Convective Tg Compari sons 
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